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1.0  INTRODUCTION 


The  primary  objective  of  this  research  effort  is  to  derive 
a  light,  maneuverable  aircraft-autopilot  model  as  one  extreme  of 
aircraft  type  for  final  approach  simulation  studies.  This  model 
is  to  serve  as  a  high  priority  vehicle  in  two  FAA  projects: 
developing  requirements  for  a  Scanning  Beam  Microwave  Instrument 
Landing  System  and  developing  an  all-encompassing  generalized 
set  of  equations  of  an  aircraft  during  approach  and  landing  for 
all-weather  landing  system  studies  on  the  NAFEC  hybrid  computa¬ 
tion  facility. 

In  general,  no  such  model  existed  for  a  common  and  repre¬ 
sentative  light  aircraft  at  the  beginning  of  this  effort  because 
light  aircraft  are  not  generally  designed  by  analysis  and  simula¬ 
tion.  Similarly,  light  aircraft  generally  do  not  have  automatic 
landing  systems  as  standard  equipment. 

The  definition  and  identification  of  a  light  maneuverable 
aircraft  is  treated  in  Reference  (1).  The  Piper  PA-30  "Twin 
Comanche"  is  selected  as  the  light  maneuverable  aircraft  for 
modeling  primarily  because  of  the  availability  of  data  from  wind 
tunnel  and  actual  flight  tests  and  because  of  an  existing, 
partially  useful,  simulation  model  at  the  NASA  Edwards  Flight 
Research  Center.  The  final  aircraft  model  is  derived  from  this 
simulation  model  and  NASA  TN  D4983.  The  flight  condition  is 
based  upon  a  high  wind  environment  (headwind  and  sidewind  of 
approximately  24  feet  per  second) .  This  condition  was  selected 
to  represent  an  extreme  case  for  the  simulation  studies. 

The  autopilot  description  is  based  on  a  report  (Reference  6) 
prepared  for  TSC  by  Dr.  Kohlman  of  the  University  of  Kansas. 

The  final  configuration  (i.e.,  gain  values,  gain  scheduling  and 
logic)  was  determined  by  simulation  at  TSC. 


2.0  VEHICLE  DESCRIPTION  4 

The  Piper  PA- 30  is  a  light  twin  engine  low-wing  monoplane. 

Figure  1  gives  the  principal  dimensions.  The  airplane  has  a 
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wing  span  of  35.98  ft.,  a  wing  area  of  178  ft.  ,  an  aspect  ratio 
of  7.3,  and  a  mean  aerodynamic  cord  of  5  feet  based  on  projec¬ 
tion  of  the  outboard  leading  edge  of  the  wing  through  the  fuse¬ 
lage.  The  wing  airfoil  section  is  a  modified  NACA642&215  air¬ 
foil  with  the  trailing-edge  cusp  faired  out.  The  wing  has  5°  of 
dihedral  with  no  twist  and  is  at  2°  positive  incidence  with 
respect  to  the  fuselage  reference  line.  The  airplane  has  the 
standard  three-control  system.  The  horizontal  tail  is  of  the 
all-movable  type  with  a  control  deflection  range  of  4°  to  -14°. 
The  tail  has  a  trailing-edge  tab  which  moves  in  the  same  direc¬ 
tion  as  the  tail  with  a  deflection  ratio  (tab  deflection  to  tail 
deflection)  of  1.5.  The  control  deflection  range  on  each  aileron 
is  from  14°  to  -18°.  The  rudder  control  deflection  range  is 
+27°. 

2.1  Vehicle  Model 

A  final  approach  model  is  presented  for  the  Piper  PA-30 
aircraft  based  on  data  available  from  the  NASA  Edwards  Flight 
Research  Center  simulation  model  and  NASA  TND  4983.  The  model 
c-nsists  of  rigid  body,  six-degree-of-freeaom  aircraft  equations 
of  motion  which  are  basically  linear  perturbation  equations  in 
the  stability  axis  system  (some  of  the  nonlinear  cross  coupling 
terms  have  been  included) . 

2 . 2  Axis  Systems 

The  stability  axis  frame  (s)  is  depicted  in  Figure  2.  The 
definition  of  airplane  angles  and  sign  convention  is  described 
in  Figure  3.  The  stability  axis  is  fixed  to  the  aircraft  and 
rotates  and  translates  with  the  aircraft.  Its  origin  is  the 
center  of  mass  of  the  aircraft.  The  X-axis  is  in  the  direction 
of  motion  of  the  airplane  in  a  reference  condition  of  steady 
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All  dimensions  are  in  feet  (meters) 

a 

Figure  1.  Piper  Dimensional  Data 
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symmetric  flight.  The  Y-axis  is  normal  to  the  aircraft's  plane 
of  symmetry  (positive  to  the  right) ,  and  the  Z-axis  is  in  the 
plane  of  symmetry  (positive  downward)  and  orthogonal  to  the  X- 
and  Z-axes.  The  aerodynamic  stability  derivatives  are  all  re¬ 
ferenced  to  this  axis  system. 

A  fuselage  referenced  body  coordinate  frame  (B) ,  is  defined 
for  determining  the  angular  rates  and  velocities  in  body  axis, 
and  the  aircraft  normal  acceleration.  This  axis  system  is 
similar  to  the  stability  axis  system  except  that  the  X-axis  is 
directed  along  the  fuselage.  The  angle  ar  relates  the  two  axis 
systems.  A  local  vertical  coordinate  frame  (E)  is  defined  for 
determining  the  velocities  of  the  aircraft  with  respect  to  the 
air  mass.  This  coordinate  frame  has  its  origin  at  the  center  of 
mass  of  the  aircraft  with  the  X-axis  pointing  North,  and  Y-axis 
pointing  East,  and  the  Z-axis  pointing  down.  The  velocities 
from  this  coordinate  frame  can  be  converted  to  ground  velocity 
by  adding  the  various  components  of  the  steady  wind. 

2.3  Assumptions  in  Using  Aircraft  Equations 

The  derivation  of  the  aircraft  equations  involved  the 
following  assumptions : 

1.  Aircraft  mass  is  constant. 

2.  The  earth  can  be  considered  an  inertial  frame. 

“  The  aircraft  is  a  rigid  body. 

4.  The  aircraft  is  symmetrical  about  its  X-Z  plane. 

5.  The  aircraft  is  initially  in  equilibrium  flight  with 
no  linear  or  angular  accelerations,  no  angular  rates, 
and  no  initial  roll  angle  or  lateral  velocity. 

6.  Small  disturbance  (perturbation)  theory  is  used. 
Motions  and  forces  are  referred  to  the  equilibrium 
flight  condition. 

Hinge  moments  are  insignificant. 


7. 


2 . 4  Equations  of  Motion 
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Pitching  Moment  Equation 
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Body  Axis  Rates  and  Velocities 


qB  =  q 

PB  =  p  cos  («r)  -  r  sin(cir) 

=  r  cos  (ar)  +  P  sin  (ar) 

UB  =  {Uo  +  U)  COS  (ar)  "  Uoa  sin  (ar) 
wB  =  U_a  cos(ar)  +  (UQ  +  u)  sin  (ar) 


Accelerations 

Vertical  acceleration  at  center  of  gravity 
w  =  UQa  cos(ar)  +  u  sin(ar) 

Vertical  acceleration  at  accelerometer 
waccel  =  WB  ”  xaccel  q 
Measured  Normal  Acceleration 
an  =  ”acoel  +  PBVB  '  qBUB 


Euler  Angle  Rates 
6  =  q  cos  $  -  r  sin  $ 

4  =  p  +  (q  sin  <t>  +  r  cos  $)  tan  9 
y  =  (q  sin  $  +  r  cos  4>)  sec  0 
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Translation  Equations  With  Respect  to  Moving  Frame 
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2.5  Vehicle  Data 


Geometry 


b  =  35.98  ft. 

c  =  5.  ft. 

S  =  178.  ft.2 

ZT  =  -.75  ft. 

aT  =  0  deg. 


m  = 


Weight 

and 

Inertias 

IxxB 

2800 

slug-f t2 

Vs  -  2801*7 

slug-ft 

IyyB  = 

1900 

slug-ft2 

Vs  =  190°- 

slug-ft 

*zz3  ~ 

4500 

slug-ft2 

I  =  4513.7 

zzs 

slug-ft 

1 xzB 

80 

slug-ft2 

I  =  -  7.9 

xzs 

slug-ft 

Trim  Flight  Condition  -  Final  Approach (High  Wind  Environment) 


=  176.  ft/sec 


.55 

.034 


(trim  elevator  setting) =  q . 4° 
(trim  flap  setting)  =0.° 


o  =  .002378  slugs/ft2 
q  =36.8  lbs/ft2 
ar  =  .0515  rad  (2.95°) 

CT  (thrust  coeff.)=  .034 

center  of  gravity  at  10%  MAC 
gear  down 
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Nor -Dimensional  Derivatives  -  Final 
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rt 
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CY6_  = 

-.00916 

2.6  Control  Wheel  and  Pedal  Characteristics 

aileron:  Gearing  constant  at  zero  control  deflection, 

0.80  radian/ft.;  wheel  deflection,  ±90°;  max. 
force  at  end  of  wheel,  ±15  lb. 

Elevator:  Gearing  constant  at  zero  control  deflection, 

0.42  radian/ft.;  wheel  throw,  4  inches  forward, 
5  inches  aft;  force  differential  40  lbs. 

(detent  position  is  4  inches  from  firewall) 

Rudder:  Gearing  constant  at  zero  control  deflection, 

0.93  radian/ft.;  pedal  deflection,  4  1/2  inches 
maximum  force  at  full  deflection,  120  lbs. 


a 


3. 


ENGINE  THRUST  DATA 


Two  160  h.p.  Lycoming  IO-320-B  four-cylinder  air-cooled  en¬ 
gines  power  the  PA-30.  Figure  4  presents  the  maximum  available 
manifold  pressure  from  each  engine  as  a  function  of  standard  tem¬ 
perature  altitude.  This  curve  should  be  mechanized  so  that  '-.he 
pilot's  indicator  shows  the  maximum  value  as  a  function  of  alti¬ 
tude  at  full  throttle. 

Figure  5  presents  the  power  available  per  engine  as  a  func¬ 
tion  of  r.p.m.,  MAP,  and  altitude.  This  curve  should  be  mechan¬ 
ized  so  that  power  cam  be  obtained  from  the  r.p.m.,  MAP,  and  alti¬ 
tude. 


Once  the  power  is  obtained,  thrust  can  be  calculated  from  the 
following  equation: 


T  (lbs.)  = 


325  Np  BHP 

(knots) 


where  N  (power  efficiency)  =  0.74 


Alternatively,  a  simplified  algebraic  expression  for  BHP  has 
been  derived  from  Figure  5  for  the  final  approach  configuration 
and  can  be  used  in  place  of  Figure  5.  The  expression  is  given  by 

BHP  =  .0024  (h)  +  .0028  (r.p.m. )  (  29 . 2- .  000989  (h)}j  -8.oJ 

where  h  =  altitude  in  feet 
ot  =  1  =  full  throttle 


An  engine  lag  of  about  0.1  second  can  be  used  in  completing 
the  engine  transfer  function  (i.e.  throttle-* thrust)  . 
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Figure  4.  Maximum  Available  Manifold  Pressure 
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AUTOPILOT 


This  section  describes  a  representative  PA-30  three-axis 
autopilot  with  stabilization  system  and  with  an  automatic  ILS  ap¬ 
proach  mode  available  down  to  decision  heiqht. 

Block  diagrams  of  the  longitudinal  and  lateral  axis  are  shown 
in  Figures  6  and  7  respectively.  Damping  is  provided  about  all 
three  axes,  and  inputs  to  the  system  come  from  a  vertical  gyro 
(6,  $),  yaw  rate  gyro  (r) ,  altitude  sensor,  and  navigation  receiver. 
Limiters  are  installed  in  the  lateral  channel  to  prevent  exces¬ 
sive  roll  angles  in  response  to  large  error  signals.  Gain  schedul¬ 
ing  is  provided  in  the  longitudinal  channel  to  desensitize  the 
system  to  glidepath  errors  as  the  runway  threshold  is  approached. 
Automatic  ILS  approaches  are  possible  down  to  decision  height. 

The  pilot  must  take  over  manually  at  that  point  because  there  is 
no  automatic  flare  capability. 

The  auto- trim  system,  which  actuates  the  elevator  trim  tab  to 
unload  the  elevator  servo,  is  not  included  in  the  block  diagrams 
since  it  has  no  measurable  influence  on  the  performance  or  dynamic 
response  of  the  airplane-autopilot  system.  It  is  a  time-delayed, 
slow,  integrating  actuator  and  as  such  does  not  respond  to  tran¬ 
sients  . 

The  output  of  the  autopilot  system  is  expressed  in  control 
surface  deflections.  Because  of  the  dependence  of  deflections  on 
dynamic  pressure,  the  final  gain  of  1/q  represents  the  aerodynamic 
gain  of  the  control  surfaces. 

The  glideslope  system  is  engaged  by  approaching  the  glidepath 
at  constant  altitude  (usually  in  the  altitude  hold  mode)  on  the 
localizer  beam.  Initially  the  glideslope  error  signal  will  be  a 
strong  nose  up  command  but  switch  B  will  be  open.  The  "glideslope 
engage  logic"  circuit  continuously  monitors  the  signal  at  point  B, 
which  will  gradually  decrease  in  strength  as  the  glidepath  is  ap¬ 
proached.  When  the  signal  at  point  B  reaches  a  0.044  rad.  (2.5°) 
nose  up  command,  switch  B  is  closed,  and  the  altitude  hold  mode  is 
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Figure  6.  Longitudinal  Control  System 
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disengaged.  To  prevent  a  sudden  pitchup  at  this  point  a  0.044 
rad-  nose  down  command  is  biased  into  the  attitude  circuit,  to 
cancel  the  input  free*  switch  B.  Thus,  after  engagement,  the  air¬ 
plane  will  continue  in  level  flight.  As  it  moves  closer  to  the 
glidepath,  the  signal  at  B  will  continue  to  decrease,  and  the  bias 
signal  will  gradually  pitch  the  nose  down.  The  gains  and  logic 
switches  are  set  so  that  the  airplane  gradually  approaches  the 
correct  glides lope  with  minimal  overshoot .  Speed  is  controlled 
manually  at  all  times  in  the  representative  autopilot  system  since 
auto-throttles  are  not  typically  provided.  However,  a  preliminary 
auto-throttle  was  derived  at  TSC  by  simulation  and  can  be  used 
for  completely  automatic  approach  studies.  The  auto- throttle  sys¬ 
tem  is  shown  in  Figure  8. 


Figure  8.  Auto- throttle 


5.  INSTRUMENT  PANEL5 

The  instrument  panel  layout  is  shown  in  Figure  9.  It  is  de¬ 
signed  to  accommodate  the  customary  advanced  flight  instruments 
on  the  left  side  in  front  of  the  pilot  and  engine  instruments  on 
the  right  side. 
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